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A model proton wave function, constructed using Poincare invariance, and 
constrained by recent electromagnetic form factor data, is used to study the 
shape of the proton. Spin-dependent quark densities are defined as matrix 
elements of density operators in proton states of definite spin-polarization, and 
shown to have an infinite variety of non-spherical shapes. For high momentum 
quarks with spin parallel to that of the proton, the shape resembles that of a 
peanut, but for quarks with anti-parallel spin the shape is that of a bagel. 



Recent data [1,2] showing that the ratio of the proton's electric and magnetic form factor 
Ge/Gmi falls with increasing momentum transfer for 1< < 6 GeV^ have created 
considerable attention. This behavior or the equivalent statement that the ratio of Pauli to 
Dirac form factors, Q F2{Q'^) / Fi{Q'^) is a approximately constant, indicates that quarks in the 
proton carry non-zero orbital angular momentum [3-6]. 

This note is concerned with the shape of the proton. Several subtleties and difficulties 
enter in determining the shape. As a particle of spin 1/2, the proton can have no quadrupole 
moment, according to the Wigner-Eckhart theorem. Thus any non-spherical manifestation 
must reside in fluctuating components of the quantum wave function. Another difficulty 
occurs in using experiments in which momentum is transferred to the proton. The proton's 
final state carries different total momentum than that of the initial state, and the effects 
of boosts cause the initial and final state wave functions to differ. For example, one often 
thinks of a particle in relativistic motion having a pancake shape because of the effects of 
Lorentz contraction. Such an effect is caused by external influence and is therefore not a 
manifestation of the intrinsic shape of the proton. Another possible problem occurs through 
the use of light cone coordinates, which involves a separation of coordinates into longitudinal 
and transverse that complicates a simple interpretation. Our technique is to use a model of 
the proton wave function [7], [5], [8], constructed with Poincare invariance and parameters 
constrained by data, to compute specific matrix elements using only the proton's rest frame 
wave function and ordinary coordinates. 

The basic experimental observables concerning us are the form factors: 



with Fi obtained from a non-spin flip matrix element, and F2 from a spin- flip term. In 
evaluating the right-hand-side of Eq. (1) we take J+ to be that of free quarks [9], 7"^ times the 
quark charge. We use a frame in which q'^ = {p' — pY, = — g^, g"*" = 0, q = Qe^ = — 2p. 

Our three-quark wave function is constructed using symmetries [11,12]. It is antisymmetric, 
expressed in terms of relative momentum variables, an eigenstate of the spin-operator defined 
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by the Pauli-Lubanski vector, is rotationally invariant and reduces to non-relativistic SU(6) 
wave function in the hmit that the quark mass goes to infinity. The wave function is given by 



^^(Pi) = $(Mo)M(pi)ti(p2)M(p3)V's(Pl,P2,P3), Pi = PiSi,Ti (2) 

where ips is a spin-isospin color amphtude factor [12], the Pi are expressed in terms of relative 
coordinates (with — psj. = K), the u{pi) are Dirac spinors and $ is a spatial wave function. 
The arguments of the spatial wave function are taken as the mass-squared operator for a 
non-interacting system [13]: 

Ma = 2(k2 + m^) Mq = (K^ + M|) + (K^ + m^) (3) 

where k is the relative momentum between quarks labeled 1 and 2 and — K is the momentum 
of the third quark [14]. The spatial wave function is of the Schlumpf [13] form: 

HMo) = ^M2^^2^, ' /5 = 0-607 GeV, 7 = 3.5, m = 0.267 GeV. (4) 

The value of 7 is chosen that Q^Gm(Q^) is approximately constant for > A GcV^.The values 
of j3 and m are determined by the charge radius and magnetic moment of the proton. Slightly 
different values of /3, 7 and m are used when the effects of the pion cloud are incorporated [8] . 
Using these newer values would cause very small differences here. 

The calculation of electromagnetic form factor is completely defined once this wave func- 
tion and operator are specified. Although the evaluation was presented long ago [7] and 
explained recently [5], it is worthwhile to briefly explain how the constant nature of the ratio 
QF2/ Fi emerges from the relativistic nature of the calculation. The wave function Eq. (2) is 
completely anti-symmetric, so we may take J"*" to act only on the third quark which absorbs 
the momentum of the virtual photon. The average charge of the third quark in the mixed- 
symmetric component of Eq. (2) vanishes, so the only component of the wave function that 
enters in the calculation of electromagnetic form factors is the mixed anti-symmetric compo- 
nent in which the first two quarks have a vanishing total angular momentum. Then the spin 
of the proton s is governed by the third quark. The relevant Dirac spinor is: 

«(p3 = K..) = ^^^f(^W-»)WV (5) 



yjE{K) + m V 



with E(K) = [K"^ -\- m?Y^'^. The total angular momentum of the proton is denoted by s and 
the lower component contains a term a ■ K that allows the quark to have a spin opposite to 
that of the proton's total angular momentum. The vector K reveals the presence of the quark 
orbital angular momentum: the struck quark may carry a spin that is opposite to that of the 
proton. Consequently helicity is not conserved [3,4,15]. 

Suppose Q is very large compared to j3 and the third quark changes its momentum from 
K to K'. Then the form factor depends on the matrix element of u{K' , s')'^^u{K, s). If we 
keep only the largest terms, those proportional to Q, we find 
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u{K',s')^+u{K,s) ~ {s'\Q + ^Qiay\s). (6) 

Mq 

The spin-flip term proportional to iay arises from the lower component of the spinor of Eq. (5), 
and this term has the same large factor Q as the non-spin-flip term. Thus Eq. (1) tells us that 
QF2 and Fi have the same dependence on Q, so their ratio is constant. The construction of an 
eigenstate of spin, mandates the use of Dirac spinors. The lower components of these spinors 
carry the quark orbital angular momentum responsible for the constant nature of QF2/F1. 

Our aim is to interpret these features of the wave function in terms of the shapes of 
the proton. The technique of deflning spin-dependent density operators is introduced by 
presenting two simple examples. Consider a non-relativistic nuclear wave function consisting 
of a proton outside a 0+ inert core. The proton is bound by a combination of a central and 
spin-orbit potential, and the single particle wave function is an eigenfunction of total angular 
momentum. Consider the case (/, js) = (1, l/2s), with one unit of orbital angular momentum. 
Then the proton wave function is a two-component Pauli spinor: {rp\ipi^i/2s) = R{rp)a ■ rp\s), 
where is the proton's position, and |s) is a Pauli spinor. The charge density of this system is 
the expectation value of the density operator 5{r — rp): with p(r) = (V'i,i/2s|^(r — rp)|'0i,i/2s) = 
R'^{r), and is spherically symmetric. However, consider instead the case that we require the 
proton at a position r to have a spin in a direction deflned by a unit vector n. Then we deflne 
a spin-dependent density: 

p(r, n) = (V'i,i/2s|(^(r - r^) ^^^^ '^^ lV'1,1/2.) 

= ^^(s|l + 2<r-f n-r-(T-n|s). (7) 

An interesting special case is to take n parallel or anti-parallel to the direction of the proton 
angular momentum s. The direction of this vector deflnes an axis (the "z-axis"), and the 
direction of vectors can be represented in terms of this axis: s ■ r — cos9. With this nota- 
tion p(r,n = s) = i?^(r) cos^^, p(r, n = — s) = i?^(r) sin^^ and the non-spherical shape is 
exhibited clearly. Note that the average of these two cases would give a spherical shape (as 
would an average over the direction of the total angular momentum), but the ability to deflne 
a spin-dependent density allows the presence of the orbital angular momentum to be revealed 
in the shape of the computed density. 

Another useful example is that of the Dirac electron wave function of the Hydrogen atom. 
This wave function is a four-component spinor given by \ipe) with 

(relV'e) = Nr2 exp{-mea Ve) ^q,/2^%^|s) ) ' 

where a is the fine structure constant and 7 = Vl — a^- We compute the expectation value 
of the spin-dependent density operator, in terms of Dirac matrices: (5(r — re))(l + 7°7-n75)/2, 
so that with n = s, p(r, n = s) = r"^''' exp (— 2meQ;re) [1 + q;^/4cos^ ^] . We see that the 
Hydrogen atom is deformed! The angular dependence is ~ 1 -|- 10~^ cos^ 9 so that the shape 
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is almost spherical. But the principle is clear: relativity as manifest by lower components of 
a Dirac wave function implies a deformed shape, if the matrix element is computed in a state 
of fixed total angular momentum. For n = — s, p(r, n = — s) = o? sin^ ^/4. 




With these examples in hand wc turn to the proton. Its wave function is specified in 
momentum space, so we define a momentum-space, spin- dependent, charge-density operator: 

P(K, n) = I (/ + 7 • 1175)^(0), (9) 

where Q/e\s the quark charge operator in units of the proton charge. The quark field operators 
are evaluated at equal time. In QCD, one would need to add a factor between the field 
operators to ensure gauge invariance. With the present model, any gluons are contained 
within constituent quarks, so the factor becomes unity. For our model, it is convenient to use 
first-quantized notation so that: 

p(K, n) = ^ {Qi/e)6{¥. - K,) (l + (7% • n75),) /2. (10) 

i=l,3 

We may compute probabilities for a quark to have a momentum K = [K, 6, (p) and spin 
direction n, for a spin-polarized proton polarized in the s direction. We find 

p(K, n) = (*,|p(K, n)|*,) = p{K)^{l + n • s + j{K){l - n • s + 2K • nK • s)) (11) 

with 

p{K) ^ J dW{k,K){E{K)+m), 7(i^) = |||y^- (12) 

4 



Some special cases of Eq. (11) are interesting. Suppose the quark spin is parallel to the 
proton spin, n = s, then p{K,n = s) = p{K) (1 + '-/{K) cos^ 6). For small K the shape is 
nearly spherical, but for large K the cos^ 6 term becomes prominent. On the other hand, 
the quark spin could be anti-parallel to the proton spin, n = — s. Then we find: p{K, n — 
— s) = p{K)'y{K) sin^ 9, and the shape is that of a torus. We may also take the quark spin 
perpendicular to the proton spin n ■ s = 0, so that p(K, n -8 = 0) = p{K){l + ■y{K))/2 + 
7(X) sin6'cos6'(cos0na; + sin^n^^), to display the dependence on the azimuthal angle. In each 
case, the non-spherical nature arises from the term proportional to 'j{K) caused by the lower 
components of the Dirac spinor (5). 




FIG. 2. Shapes of the proton, s is in the vertical direction z. Left column quark spin parallel to 
nucleon spin n = s. Right column : quark spin anti-parallel to nucleon spin n = — s. The value of 
K increases from to 1 to 4 GeV/c. 

We turn to the numerical evaluation and display of these shapes. The shape for a given 
value of K is determined by the ratio 'y{K) which reaches a value of 0.6 for K = 1 GeV/c. 
This implies considerable non-sphericity. The probability that a given value of K is reached 
is determined by the function K'^p{K), displayed in Fig. 1. The most likely value of K 
corresponds to 7(-ft') = 0.16, and ■y{K) < 1. 

The shapes for the cases of quark spin parallel and anti-parallel to the polarization direction 
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of the proton s are displayed in Fig. 2. As the value of K increases from to 4 GeV/c the 
shape varies from that of a sphere to that of a peanut, if n || s. The torus or bagel shape is 
obtained if — n || s. Taking n ± s leads to some very unusual shapes shown in Fig. 3. 




FIG. 3. n • s = 0. Left column, n = x (out of page), central: n = y, right n = (x + y)/^/2. The 
momentum K increases from 1 to 4 GeV/c. 

It is worthwhile to obtain a coordinate space look at the shape of the proton. One can 
determine the probability for the separation of a quark from the center of mass to be r because 
r is canonically conjugate to K. We define a coordinate-space, spin- dependent, charge-density 
operator that acts on the quarks in the cm frame: 

/)(r, n) = 5: {Q^/e)5{v - r,)(l + (7% " n75).)/2 (13) 

i=l,3 

so that we may compute probabilities for a quark to have a position r = (r, 6, (p) and spin 
direction n. We find 

p(r,n) = (*,|p(r,n)|^,) = y kx\k, r)^{l + -f^^-f 11^5) x{k,r), (14) 

with 

and Fk{r)= J d'K^k,K)iE{K)+My/^e^^-^,Gk{r) = f J d^K^k, K)j^^^^^^. We find 

p(r, n) = pc7(r)i (1 + n ■ s) + pL(r)i (1 + 2f • sn ■ r - n ■ s) , (16) 

where pu{r) = / d'^kF^{r), phij") = / d^kG\{r). The pattern is similar to the one in mo- 
mentum space, with p(r, n = s) = puij") + piir) cos^ 9, p(r,n = — s) = pi^(r) sin^ 6^, and 
p(r, n = ^) = |pt/(r) + |pL(r) (l + 75 cos 6 sin 6'(cos + sin <; 
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FIG. 4. Coordinate space densities 

The ratio Pl/ Pu controls the relativistic effects. This can be much larger (Fig. 4) than 
the factor 7(i^) (limited by unity) controlling the momentum-space shapes, so more extreme 
shapes are possible. The most likely value of PlI Pu is about 0.25, but there is no limit. The 
case with pu/ Pu=^ is shown in in Fig. 5. A pretzel form is obtained if n is out of the page. 




FIG. 5. Shape of the proton coordinate space. Left: n = s, right: n points out of the page 

The shape of the proton may be defined in terms of matrix elements of spin-dependent 
density operators Eqs. (10,13) taken for protons in any fixed polarization state. Relativity 
mandates the use of Dirac spinors to describe the quarks. These components, embodied in 
Eq. (6), lead to a constant ratio of QF2/F1 in accord with observation, and also to shapes of 
an infinite variety, Eqs. (11,16), Figs. 2,3,5. 

We next consider experiments capable of measuring the matrix element /3(K, n), Eq. (9), 
for real nucleons |A^). Observe that / d?K p{K.,n) is a local operator. Its matrix element is a 
linear combination of the charge, integrals of spin- dependent structure functions Ag, and qa, 
it can be determined from previous measurements. We find 

/ d'K{N\p{¥., n = ±s)|iV) = hN\m-{l' ± 7V)V^(0)|iV) 
J 2 e 
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= ^ (l ± {^{Au + Ad + As) + ^QA^ = 0.5 ± 0.34, (17) 

in which numerical values of Aq are taken from Ref. [16]. The model we use gives 0.5 ± 0.37 
for the above quantity, indicating that the results shown here may not be unrealistic. 

The task of determining p(K, n) as a function of K remains open. However, the specific 
relativistic effects of orbital angular momentum and the related spin-flip effects responsible 
for the non-spherical shapes can be expected to influence many measurable quantities. These 
include: spin-dependent structure functions gi,g2, the NA electromagnetic transition form 
factor, and spin-dependent cross sections in high momentum transfer scattering of polarized 
protons. These topics will be addressed in other publications. 
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